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We explore the zero-temperature phase diagram of bosons interacting via Feshbach resonant 
pairing interactions in one dimension. Using DMRG (Density Matrix Renormalization Group) 
and field theory techniques we characterize the phases and quantum phase transitions in this low- 
dimensional setting. We provide a broad range of evidence in support of an Ising quantum phase 
transition separating distinct paired superfluids, including results for the energy gaps, correlation 
functions and entanglement entropy. In particular, we show that the Ising correlation length, order 
parameter and critical properties are directly accessible from a ratio of the atomic and molecular 
two-point functions. We further demonstrate that both the zero-momentum occupation numbers 
and the visibility are in accordance with the absence of a purely atomic superfluid phase. We 
comment on the connection to recent studies of boson pairing in a generalized classical XY model. 
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I. INTRODUCTION 

The rapid progress in manipulating ultra cold atomic 
gases has led to new approaches to strongly interacting 
quantum systems. This includes the properties of highly 
correlated states of matter, such as Bose-Einstein con- 
densates (BECs) [U [2], Mott insulators [3], and super- 
solids It also allows access to the phase transitions 
and crossovers between these fascinating phases. In the 
last few years, the BEC-BCS crossover [5J{5] between a 
molecular BEC and a Bardeen-Cooper-Schrieffer (BCS) 
pairing state, has stimulated a wealth of experimental 
activity using fermionic atoms [SHIS]- This has been 
achieved through the use of Feshbach resonances, which 
enable one to control the strength of pairing interactions 
using a magnetic field. This has not only opened the 
door to central problems in condensed matter physics, 
but also offers insights into the quantum chemistry of 
molecule formation and chemical reactions jl7) . 

In tandem with these advances, Feshbach resonances 
and molecule formation have also been studied in bosonic 
systems [18 26^. Recent experiments have been per- 
formed both in optical traps [?7l - B5] and in optical lat- 
tices [3D]; for a review see Ref. [3T]. On the theoretical 
side, the BEC-BCS "crossover" problem for bosons has 
also been investigated in the continuum limit [32 - 34 and 
on the lattice [531110] . The problem differs markedly from 
the fermionic case since the carriers themselves may Bose 
condense. This leads to the possibility of an Ising quan- 
tum phase transition occurring between distinct paired 
superfluids |5^H51] . The phases are distinguished by the 
presence or absence of carrier condensation, and the as- 
sociated quantum phase transition involves discrete Z2 
symmetry breaking. Closely related phases and quan- 
tum phase transitions have also been observed in multi- 



component fermion systems |41H44j . and in the attrac- 
tive Bose-Hubbard model with a restricted three parti- 
cle Hilbert space [45l - f5T] . There are also magnetic ana- 
logues in quantum spin chains [S^]. More recently, the 
phenomenon of boson pairing has also been explored in 
the context of a generalized classical XY model with two 
competing harmonics in the periodic interactions [53 - 56] . 
This has led to the prediction of a novel phase diagram 
with unusual criticality. 

Motivated by the possibility of stabilizing pairing 
phases of bosons [57H5^ in cold atomic gases, we recently 
investigated the bosonic Feshbach resonance problem in 
one dimension (ID) |63j. We employed large scale DMRG 
(Density Matrix Renormalization Group) |M] and field 
theory techniques [63) in order to incorporate the effects 
of enhanced quantum fluctuations in ID. Amongst our 
findings, we presented compelling evidence for an Ising 
quantum phase transition separating distinct superflu- 
ids [5S]. The aim of the present manuscript is to shed 
further light on this novel transition, and to provide a 
thorough discussion of the superfiuid phases in this ID 
setting. In particular, we describe a variety of meth- 
ods to extract the Ising characteristics from the gapless 
superfluid background. We also provide a quantitative 
finite-size scaling analysis of the zero-momentum occu- 
pation numbers and the visibility. Our results are con- 
sistent with the absence of a purely atomic superfluid 
phase with non-condensed molecules, in contrast to the 
earlier suggestions of Ref. [371 ISB] • 

The layout of this paper is as follows. In Sec. |ll| 
we present the Hamiltonian under investigation and in 
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Sec. Ill we discuss the phase diagram, 
describe the associated field theory and gather our pre- 
dictions for a variety of local expectation values and cor- 
relation functions. We use these results to characterize 
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the different pliases and to establisfi a detailed compar- 
ison with DMRG. In Sec. |V] we provide a quantitative 
account of the finite-size scaling of the zero-momentum 
occupation numbers and the visibility. We contrast our 
results with those of Refs. [371 HI]- In Sec. VI we dis- 
cuss the behavior of the entanglement entropy and the 
emergence of Ising criticality at the transition between 
the distinct paired superfluids. We also discuss the be- 
havior at the superfluid-Mott insulator transitions. In 
Sec. |VII| we describe the Ising scaling regime, and dis- 
cuss a variety of ways to extract the principal Ising char- 
acteristics. This includes the Ising order parameter and 
the correlation length using a finite-size scaling analysis 
of the atomic and molecular correlation functions. We 
also discuss the utility of a suitable ratio of the atomic 
and molecular two-point functions for analyzing the Ising 



quantum phase transition. We conclude in Sec. VIII and 
provide further directions for research. 



II. MODEL 

We consider the Hamiltonian 



'i a 



(1) 



describing bosons, bia, hopping on a ID lattice with sites 
i, where a = a,m labels atoms and molecules 35-40J. 
Here, are on-site potentials, ta are nearest neighbor 
hopping parameters, and Uaa' are interactions. We as- 
sume that molecule formation is described by the s-wave 
Feshbach resonance term. 



m\aiai + H.c), 



(2) 



where we denote mrii = and = bia] for re- 

cent work on the p-wave problem see Refs. [BSJ 1^ . 
This conversion implies that the number of atoms and 
molecules are not separately conserved, but the total, 
A^T = ^ii^ia + 2ni„), is preserved. For simplicity, in 
writing Eq. ([T]) we neglect any effects of higher Bloch 
bands in optical lattices [68l - f70] . In this respect, the 
Hamiltonian ([T]) may be regarded as a lattice regular- 
ization of the continuum models studied in Refs. [5^H5i] : 
see also Refs. [^^26) . This approach is very convenient 
for numerical simulations, and enables us to investigate 
the superfluid transitions where lattice effects are ger- 
mane. It also allows us to make contact with previous 
quantum Monte Carlo (QMC) simulations [37, 38J and to 
place the problem on a firmer footing. As in the original 
works , we neglect the effects of three body losses 

and finite molecular lifetimes. 

In this manuscript we use DMRG on ID systems with 
up to L = 512 sites, where we set the lattice spacing to 



unity and adopt energy units where ta = \. We further- 
more set tm = 1/2 throughout. We work in the canon- 
ical ensemble with the total density — N^^/L = 2 
held fixed and allow up to five atoms and five molecules 
per site, corresponding to a large Hilbert space of dimen- 
sion (6 X 6)^; for a discussion of the effects of chang- 
ing the local Hilbert space dimension see Appendix \K\ 
With open (periodic) boundary conditions we retain up 
to TOp = 2400 {nip = 3000) states in the density matrix 
in order to ensure that the discarded weight is less than 
1 X 10-1° (1 X 10-S). 



III. PHASE DIAGRAM 

As we discussed in Ref. 'S3], the qualitative phase di- 
agram of the ID lattice Hamiltonia n |l I was previously 



considered using QMC simulations |37l [38| . In addition 
to delineating the Mott insulating and superfluid phase 
boundaries, this work led to intriguing predictions of su- 
perfluidity within the Mott phase, and an additional su- 
perfluid phase not present in mean field theory |35H33]. 
Although we find very good quantitative agreement with 
many of the numerical results [371 EH] , these additional 
predictions are at variance with our recent findings [63] 
which combine field theory with DMRG. This was also 
suggested by our earlier studies using hardcore bosons 
[391 |40] . It has recently been argued that the absence 
of particle conservation hindered the interpretation of 
these previous QMC simulations |71j . In this manuscript 
we will further demonstrate that the use of momentum 
space observables, including the zero-momentum occupa- 
tion numbers and the visibility, also complicated the in- 
terpretation of these earlier finite-size QMC simulations. 

In order to put the problem on a more stable platform, 
we present a section of the phase diagram in Fig. [l] with 
parameters chosen for comparison with Ref. [37j- (Note 
that our conventions differ from Ref. |37' by a factor of 
1/2 in the interaction terms so that double occupancy 
corresponds directly to Uaa' ■ Also, e^n plays the role of 
their detuning parameter, D, when — 0.) The phase 
boundaries shown in Fig.[l]correspond to the vanishing of 
the one-particle and two-particle excitation gaps, Eig = 
f^ip{L) - fiih{L) and £^29 = fJ'2p{L) - H2h{L) respectively, 
where the data are extrapolated to the thermodynamic 
limit. L — !■ 00. Here 



finp{L) = [Ea{L,NT + n)-Eo{L,NT)] /n, 
finhiL) = [Eo{L,Nt) - Eo{L,Nt - n)] /n, 



(3) 



where Eq{L,N) is the ground state energy for a system 
of size L and a total number N of atoms and molecules. 
The phase diagram in Fig. [T] consists of three distinct 
phases: a Mott insulator (MI) with gaps for both exci- 
tations, Eig ^ and E2g ^ 0, a molecular condensate 
(MC) with a one-particle gap Eig ^ and i?2g — 0, and 
a coupled atomic plus molecular condensate (AC-I-MC) 
with Eig = and E2g = 0. As we shall discuss more 
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FIG. 1. (color online), (a) Phase diagram of the ID Hamilto- 
nian ([T]) with total density pT = Nt/L — 2, showing a Mott 
insulator (MI), a molecular condensate (MC), and a coupled 
atomic plus molecular condensate (AC+MC). We use DMRG 
with up to L = 128 sites and open boundary conditions. We 
choose parameters ta = 0, Uaa/2 = Umm/^ = Uam = g = U, 
ta = 1, tm = 1/2, for comparison with Ref. [ST]- The squares 
and circles indicate the vanishing of the one-particle and two- 
particle gaps, Eig and £29, respectively, as L — ^ 00. The stars 
and crosses indicate where the molecular and atomic correla- 
tion exponents, Vm and Va reach 1/4 in the MC and AC-(-MC 
phases respectively. These values correspond to a molecular 
KT transition and an atomic KT transition respectively. The 
remaining panels (b), (c) and (d) show the variation of the 
extrapolated gaps, Eig and i?2g with L — > 00, on passing 
through the phase boundaries at the corresponding points in 
panel (a). In particular, we provide detailed evidence for an 
Ising quantum phase transition occurring between the MC 
and AC-f MC phases. 



fully below, the MC phase may be interpreted as a pair- 
ing phase of bosons in the absence of atomic condensa- 
tion. In contrast, the AC-I-MC phase has both molec- 
ular and atomic condensation. In comparison to the 
qualitative phase diagram presented in Ref. 07] , inferred 
from quantum Monte Carlo simulations on smaller sys- 
tem sizes, we find no evidence for a single-component 
atomic superfluid phase co-existing with non-condensed 
molecules. This is in accord with theoretical expecta- 
tions in higher dimensions, where atomic condensation 
is always accompanied by molecular condensation [32- 
133] provided the molecules are present; in the extreme 
limit where 6^ — ^ 00, occurring on the boundary of the 
AC-I-MC phase, the molecules are explicitly excluded by 
the chemical potential as shown in Fig. [2j The conclu- 
sions of Ref. [37] have also come under scrutiny due to 



the additional claims of superfluidity within the Mott 
phase [39l[7T]. Here, however, our main focus is on the 
character of the transition between the distinct MC and 
AC-I-MC superfluids. In the subsequent discussion we 
will begin with symmetry arguments and field theory pre- 
dictions before turning to a comparison with DMRG. 



IV. FIELD THEORY DESCRIPTION 

A heuristic way to understand the possibility of an 
Ising quantum phase transition between the distinct MC 
and AC-I-MC superfluids is via the generic number-phase 



relationships, a ~ ^/pae and m ~ ^Pm e , where pa 
and pra are the average atomic and molecular densities 
respectively. Substituting these expressions into ([T]), the 
Feshbach term ([2]) takes the form [51] 



By 



'^9 Pay 



cos(z?™ - 2'da). 



(4) 



Minimizing this interaction locks the phases of the atomic 
and molecular condensates according to the relationship 



(5) 



where for simplicity we assume g > 0. We see that the 
phases are locked, but only modulo tt, and this gives 
rise to the possibility of a discrete symmetry breaking Z2 
transition between Feshbach coupled superfluids. Denot- 
ing -dm = "i?, one may recast the number-phase relation- 
ships in the form |34j 



(6) 



where the Feshbach locking is explicitly enforced and 
^ \fPo.^^^^'^ plays the role of an Ising degree of free- 
dom. The decomposition ([6| will play a central role in 
the subsequent analysis and allows one to gain a handle 
on the correlation functions and the principal features of 
the phase diagram. 

An alternative way to understand the possibility of 
an Ising quantum phase transition between the MC and 
AC+MC phases is via the symmetry of the Hamiltonian 
([!]) under U(l) x Z2 transformations: 



m — e m, 



(7) 



where G M. Before discussing the problem in ID, where 
continuous symmetry breaking is absent, we first con- 
sider the behavior in higher dimensions j32H34j . In this 
case, the molecular condensate (MC) phase has (m) 7^ 
and (a) — 0. This only breaks the U(l) symmetry, 
and leaves the Z2 symmetry, a — > —a, unbroken. This 
corresponds to an Ising degree of freedom in the disor- 
dered phase, coexisting with molecular superfluidity. In 
contrast, the coupled atomic plus molecular condensate 
(AC+MC) phase has (m) ^ and (a) ^ 0. This breaks 
the U(l) X Z2 symmetry completely and corresponds to 
a Z2 ordered Ising degree of freedom, coexisting with 
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atomic and molecular superfluidity. Returning to the 
present ID problem, where continuous U(l) symmetry 
breaking is prohibited, the formation of expectation val- 
ues for (a) and (m) is excluded. Instead, superfluidity is 
characterized by power-law correlations, and the nature 
of the phases and quantum phase transitions in Fig. [l] 
requires further examination. 

Due to the U(l) x Z2 symmetry of the Hamiltonian 
Q, the low energy Lagrangian of the MC to AC-I-MC 
transition is given by £ = £^ -I- £0 -f [Ml [l2] , where 



(8) 



is a free bosonic field and 



c-^\dr4>Y + [d^^ct^Y ~M4>' + W (9) 



is an Ising model in the soft-spin (/i* representation. The 
coupling between the two sectors, C^^ = i4P'dr'&/'i, has a 
form similar to a Berry phase [Ml [T^ ■ A closely related 
action also arises for tunnelling between quantum wires 
[73) . In the following we will neglect the contribution C^^ 
and explore the consequences of the reduced action. Suf- 
ficiently far away from the transition this can be justified 
by a mean-field decoupling, C^^ ~ i{(j))^dr'd/2, which 
reduces the additional interaction to a total derivative 
term, which can be neglected. The simplified action is 
therefore expected to provide a good description of the 
proximate phases. Near the quantum phase transitions, 
this cannot be neglected a priori, and £,30 may change 
the behavior on very large length scales and in other re- 
gions of the phase diagram [75]. However, all of our find- 
ings are consistent with expectations based on + 
only. The parameters K^, c^, if^, c^, 77, A, are related 
to the coefficients of the Hamiltonian (Tj), but the details 
need not concern us here. In this field theory approach, 
the atoms and molecules are described by the semiclas- 
sical number-phase relations given in Eq. In the 

subsequent discussion we will explore the ramifications 
of this correspondence in ID, both for local observables 
and correlation functions. For complementary work using 
the Bethe Ansatz and bosonization see also Refs. [7411751. 
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FIG. 2. (a) Local expectation values obtained by DMRG 
for the ID Hamiltonian (|T]) with open boundaries and up to 
L = 128 sites. We use the same parameters as in Fig. [l] 
and set U = 0.7. We plot the local density of atoms (ajai) 
(circles), molecules {mlrrii) (squares), and the expectation 
value of the Feshbach conversion term {m\aiai) (triangles), 
evaluated at the system mid-point, i = L/2. All of these 
quantities are generically non-zero on both sides of the MC 
to AC-I-MC transition as indicated by the dashed line. In 
the limit of large positive (negative) detuning we have 
mainly atoms (molecules) with a density determined by the 
canonical ensemble constraint pT ~ 2. (b) Corresponding 
finite-size data and linear extrapolation as a function of 1/L 
for em = —3 (open) and e,n = —5 (filled). 



A. Local Expectation Values 

An immediate consequence of the decomposition Q is 
that the densities of atoms and molecules 

{m^x)m{x)} ^ p,n, {a\x)a{x)} ^ {(jY), (10) 

are generically non-zero in both the MC and AC-I-MC 
phases. This is supported by our DMRG results as 
shown in Fig. [2]ja). These are extrapolated from the 
finite-size data to L — >■ 00, as indicated in Fig. [2][b). 
In the limit of large positive detuning with — >■ cxd 
we have mainly atoms as one would naively expect and 
{a^ {x)a{x)) ^ 2. Likewise, in the limit of large negative 



detuning, corresponding to e„i — >■ — cx3, we have mainly 
molecules and (m^(x)m(a;)) ~ 1. These limiting densi- 
ties are consistent with working in the canonical ensemble 
with pt = X^il^-ia + 2nim)/i = 2 held fixed. 

In addition to these local densities the local expecta- 
tion value of the Feshbach conversion term, 

{mHx)a{x)a{x)) ^ ^{4>'') ^ 0, (11) 

is non-zero. It exhibits true long range order, even in 
this low-dimensional setting. This is a consequence of 
the relevance of the Feshbach term in the renormalization 
group sense. Our numerical results in Fig. |2] show that 
this quantity is indeed finite. In particular, this confirms 
the locking of the phases of the atomic and molecular 
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condensates (modulo tt) on both sides of the transition. 
However, due to the symmetry under a — —a, the expec- 
tation value ( [TT] ) is naively insensitive to the Ising transi- 
tion itself, as may be seen in Fig. [2] Further insight into 
this quantum phase transition and the proximate phases 
is more readily obtained from correlation functions. We 
will explore this in more detail below. 



B. Green's Functions and Pairing Correlations 

The nature of the MC and AC-I-MC phases shows up 
most clearly in the atomic and molecular Green's func- 
tions, {a^x)a{0)) and {m^x)m{0)), and the pairing cor- 
relations {a' {x)a^ {x)a{0)a{0)) . Their spatial dependence 
is dictated by the correlations of the underlying Ising 
model in Eq. ([9]), and we address each phase in tiun. 

1. Z2 Disordered MC Phase 

As follows from the decomposition (|6]), the molecular 
Green's function 



(12) 



decays as a power-law, where the correlation exponent, 
i^m = ^/{2ttK§) varies throughout the phase diagram, 
and ao is a short-distance cutoff. In contrast, in the Z2 
disordered MC phase, the atomic Green's function 

{a\x)a{0)) (X (0(:E)e-^^0(O)e^^) ^ (^)^Ko(x/C), 

decays exponentially, where ^ is the Ising correlation 
length. Here we use the hard-spin fermionic represen- 
tation of the Ising model to write [75] 



Wx)m) - Ko(x/0, 



(13) 



where Kg is a modified Bessel function. On the other 
hand, pairs of atoms exhibit power-law correlations 



a\x)a^x)aiO)a{0)) 



X 



(14) 



where the exponent i^i, ~ for these atomic bilinears 



coincides with the molecular exponent in Eq. (12 1. That 



is to say, the MC phase is a pairing phase of bosons with- 
out power-law atomic condensation [57H62] . 

In order to explore these field theory predictions in 
more detail we perform DMRG on the ID Hamiltonian 
([ij). The predicted behavior is well supported by our 
simulations in Fig. [3[ The molecules and atomic bilin- 
ears show power-law behavior with the same exponent, 
'^m — i^b, whilst the atomic two-point function shows ex- 
ponential decay. Our DMRG results also indicate that 
this behavior persists into the regime close to the Mott 
insulating phase boundary shown in Fig.[T] In particular, 
the molecular correlation exponent reaches the value of 
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FIG. 3. Correlation functions in the Z2 disordered MC phase 
obtained by DMRG on the ID Hamiltonian ([l] with L = 512 
and open boundaries. Here and throughout the manuscript 
we consider sites displaced around the system mid-point in 
order to minimize boundary effects. We use the same pa- 
rameters as in Fig. [ij and set U — 0.7 (open), U = 0.5 
(filled), Em ~ —4 (circles) and e,„ = —6 (triangles), (a) 
Atomic Green's function {a\aj) showing exponential decay, 
(b) Molecular Green's function {m\mj) showing power-law 
behavior, (c) Bilinears of atoms {a\a\ajaj) showing power- 
law behavior with the same exponent as the molecular Green's 
function in panel (b); see Fig. [4] This establishes the MC 
phase as a pairing phase of atoms without power-law atomic 
condensation. 



Vm — 1/4 at the MI boundary; see Fig. |4] This is consis- 
tent with a molecular Kosterlitz-Thouless (KT) [771 ES] 
transition. It is analogous to the behavior at the tips of 
the Mott lobes in the single-band Bose-Hubbard model 
[79H8T] where the Luttinger liquid parameter takes the 
value K — 1/{2'kv) — 2/7r in the normalization con- 
ventions of Eq. ([8|. The latter transition takes place 
at constant density, and is therefore compatible with our 
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FIG. 4. (color online). DMRG results for the two-point func- 
tions of (a) molecules and (b) atomic bilinears in the MC 
phase with open boundaries and up to L = 512. We extract 
the molecular and bilinear exponents, i/m and Vb, by finite-size 
scaling collapse of the data for different system sizes. In (c), 
(d) and (e) we show the resulting evolution of Vm (circles) 
and ut (stars), for vertical scans through Fig. [l] with fixed 
values of e™. The vertical dashed lines correspond to the lo- 
cation of the MC to MI transition obtained from the gap data. 
The molecular exponent reaches the value of Urn = 1/4 at the 
MC to MI transition. This corresponds to a molecular KT 
transition and is analogous to the fixed density transition at 
the tips of the Mott lobes in the single-band Bose-Hubbard 
model. The critical exponent Vb associated with the power- 
law decay of the atomic bilinears {a\a\ajaj) (stars) coincides 
with Vm- 



canonical ensemble constraint, = 2. 

We recall that in deriving the above correlation func- 
tions we have neglected the coupling term, in the 
low-energy Lagrangian so that the expressions factorize 
into independent U(l) and Z2 contributions. The good 
agreement with DMRG lends a postiori support to this 
approximation within the explored region of the phase 
diagram. 



2. Z2 Ordered AC+MC Phase 



In the Z2 ordered phase the molecular Green's function 
(mt(x)m(O)) oc (e-"'^^) e'-'^o)) ^ (^)''" , (15) 



continues to decay as a power-law. 
atomic Green's function 



In addition, the 



'■2 C 



V X 



(16) 



also decays as a power-law, where the atomic correlation 
exponent i>a = ^'m/4 is locked to the molecular exponent 
by a factor of one quarter [72]. This is a consequence 
of the Feshbach coupling which ties the phases of the 
atomic and molecular condensates together. Note that in 
writing Eq. (16), we approximate the result for the two- 



point function of the Ising order parameter at leading 
order 



where 



(0)2[l + 7r-2F(x/e)] (17) 



F(z) = z^[Kiiz)~Kliz)]~zKoiz)K,iz) + -Kliz), (18) 

and Ko(z) and Ki(z) are Bessel functions. These pre- 
dictions of power-law behavior, as given by Eqs. ( 15 ) 



and (16), are well supported by our numerical simula- 
tions as shown in Fig. [5j The locking of the atomic and 
molecular correlation exponents Va = is also ob- 

served. In addition, these robust features persist into 
the large-?/ regime where field theory arguments are no 
longer strictly valid. In particular, the atomic and molec- 
ular correlation functions remain as power-laws right up 
to the MI boundary shown in Fig. [T] We find that the 
atomic exponent i>a reaches the value of I'a = 1/4 at the 
AC+MC to MI transition; see Fig. [6j This is consistent 
with an atomic KT transition, as occurs at the tips of 
the Mott lobes in the single-band Bose-Hubbard model. 
At the same time, the molecular exponent reaches 
the value of Vm = 1 due to the aforementioned exponent 
locking; see Fig. [6j The presence of this molecular su- 
perfluid close to the MI boundary, clearly supports the 
absence of a single component atomic superfluid phase in 
this ID setting, in contrast to the findings of Ref. [37] , 
This is also compatible with mean field theory in higher 
dimensions |32H34j where atomic condensation is always 
accompanied by molecular condensation due to the struc- 
ture of the Feshbach term. Hp. We will return to this 
issue in Sec. [V] in our discussion of the corresponding 
zero-momentum occupation numbers and the visibility. 



3. Mixed Correlation Functions 

In addition to the purely atomic or molecular Green's 
functions, it is also instructive to examine the mixed cor- 
relation functions involving both atoms and molecules. 
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FIG. 5. Correlation functions in the Z2 ordered AC+MC 
phase obtained by DMRG on the ID Hamiltonian |TJ with 
L — 512 and open boundaries. We use the same parameters 
as in Fig.[l]and set U = 0.7 (open), U = 0.5 (filled), = -2 
(circles) and = —3 (squares), (a) Atomic Green's function 
(alaj) showing power-law decay, in contrast to Fig. da), (b) 
Molecular Green's function {m\mj) showing power-law be- 
havior. The exponent tracks the atomic exponent in (a) up 
to a factor of 4; see Fig. [g] (c) Bilinears of atoms l^a\a\ajaj) 
showing power-law behavior with the same exponent as the 
molecular Green's function in panel (b). This establishes the 
AC-I-MC phase as a pairing phase of atoms in the presence of 
atomic condensation. 




0.5 



1 

3.4 


3.6 


1 . 1 1 
3.8 4 


1 ' 1 ' ' 

1 


(d) e,„ = 






1 




AC+MC 1 MI - 
1 

1 1 , 


2.5 




3 



0.5 



' 1 ' + 


(e)e,„ = -2 








» * ® * 


* AC+MC 1 MI - 
1 , 1 



1.5 



2.5 



U 



FIG. 6. (color online). DMRG results for the two-point func- 
tions of (a) atoms and (b) molecules in the AC+MC phase 
with open boundaries and up to L = 512. We extract the 
atomic and molecular exponents, Va and Um, by finite-size 
scaling collapse of the data for different system sizes. In (c), 
(d) and (e) we show the resulting evolution of Aua (circles) and 
Vm (stars), for vertical scans through Fig. [l] with fixed values 
of em- The vertical dashed lines correspond to the location 
of the AC+MC to MI transition obtained from the gap data. 
The data confirm the locking of the atomic and molecular ex- 
ponents via the relation Vm = '^Va- The exponents reach the 
values of i^a = 1/4 and Um = ^ sX the MI boundary. This is 
consistent with an atomic KT transition. It is analogous to 
the behavior at the tips of the Mott lobes in the single-band 
Bose-Hubbard model. 



It follows from the decomposition ([6| that 

(mt(x)a(0)a(0)) ^ ^{^^) (^)'" , (19) 

decays as a power-law with the same exponent as the 
molecular Green's function. Once again, this reflects the 



phase locking of the atomic and molecular condensates 
due to the Feshbach term, and is present in both the MC 
and AC-I-MC phases. This behavior is in very good agree- 
ment with our DMRG simulations as shown in Fig. [T] 
In particular, the power-law exponent tracks those dis- 
played in Figs. |3|b) andjsjjb) for {m){x)m{Q)). 
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FIG. 7. DMRG results for the mixed correlation function 
— {m\ajaj) with L = 512 sites, open boundaries and U — 0.5. 
The data correspond to e^n ~ —2 (circles), = —3 (squares), 
em = —4 (up triangles), em ~ —5 (diamonds), e„i = —6 (down 
triangles) and show power-law behavior in both the AC and 
AC+MC phases. As predicted by Eq. (191, the exponents 



agree with those of the molecular Green's function in panels 
(b) of Figs. [3] and [5] 



C. Density Correlation Functions 

Having discussed the atomic and molecular Green's 
functions we now turn our attention to the correlation 
functions of the local densities. Denoting nm{x) = 
{x)m{x) and na{x) = a^{x)a{x) one obtains 



(x) ^ Pm+ lldx'd + . ■ . , 



(20) 



where pm and pa are the average molecular and atomic 
densities, and 71, 72, 73 are constants. Here we use the 
primary correspondence given in Eq. ([6]), and combine 
the exponentials by point-splitting and the short distance 
operator product expansion. The expansion (20) incor- 



porates the effects of density fluctuations and it follows 
that the density-density correlations have the same lead- 
ing dependence in both the MC and AC-f MC phases: 



{na{x)ni3{0)) ~ PaPfS 



Call 

x^ 



(21) 



where a,/3 £ a,m and Cap are non-universal constants. 
This is confirmed by our DMRG results in Figs. [8] and |9] 



V. MOMENTUM SPACE OBSERVABLES 

In the previous section we have focused directly on 
the superfluid correlation functions due to the absence 
of continuous symmetry breaking in ID. However, a use- 
ful diagnostic of superfluidity in higher dimensions is the 
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FIG. 8. DMRG results for the connected density correlation 
functions in the Z2 disordered MC phase for the parameters 
used in Fig. [T] The values of em are indicated in panel (b) . 
We use open boundaries with L — 128 and set U — 0.7. 

(a) [{n^anja) - {nia){nja)\. (b) \{nimnjm} - {nim}{njm)\- (c) 
l(n-mWjm) — {'nia){njm)\- The results are in a gree ment with 



the leading 1/x dependence predicted by Eq. (21 1 



divergence of the occupation number 



,(fc) 



1 Y: e^'^'-'H^.a^^,) 



(22) 



at zero momentum, fc = 0. This quantity was recently 
used in Ref. |37| . in conjunction with visibility data, to 
argue in favor of a single component atomic superfluid 
phase in the ID system ([T]). In view of our results in 
the previous sections, which show the presence of both 
atomic and molecular superfluidity right up to the Mott 
boundary in Fig.jl] we revisit this issue here. As shown in 
Fig. [6j the atomic and molecular correlation functions in 
the AC-I-MC phase are power-laws, {m} {x)m{0)) ^ x^"™ 
and {a^ {x)a{0)) ~ a;"'''"/'*, with locked exponents. Sub- 



stituting these asymptotic forms into Eq. ( 22 1 suggests 
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FIG. 9. DMRG results for the connected density correlation 
functions in the Z2 ordered AC+MC phase for the parameters 
used in Fig. [T] The values of are indicated in panel (b). 
We use open boundaries with L — 128 and set U — 0.7. 
(a) {{niarija) - {nia){nja}\- (b) \{nimn-im) ~ {nim}{njm)\- (c) 
\{nianjm) — {nia){njm,)\- The results are in a gree ment with 
the leading dependence predicted by Eq. (21 1. 



that the zero-momentum occupation numbers depend on 
system size according to [Ml [51] 



(23) 

where Aa^m and Sa,m are constants. In particular, since 
the molecular exponent, z/„j, only reaches unity at the 
Mott phase boundary (see Figs, [l] and [6]) both of these 
zero-momentum occupation numbers are expected to di- 
verge with increasing system size. This is supported by 
our DMRG results as shown in Fig.[TO] However, it is ev- 
ident from Eq. ( 23 1 that (0) diverges very slowly with 
increasing system size close to the MI boundary since 
— i' 1. In the absence of a detailed finite-size scal- 
ing analysis this may lead to the erroneous conclusion 
of a purely atomic superfiuid. In addition, our findings 
suggest the absence of any change in behavior in the con- 
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FIG. 10. DMRG results for the dependence of the zero- 
momentum occupation numbers na(0) and nm(0) on system 
size L within the AC-I-MC phase shown in Fig. [l] The results 
are consistent with algebraic correlations for both atoms and 
molecules with locked exponents Vj^ — Ava- The presence of 
molecular superfluidity in the lower panels confirms the ab- 
sence of an AC phase close to the MI boundary. 



vergence properties of rim(O) as L — )■ 00, which could be 
misinterpreted as a quantum phase transition to a purely 
atomic superfiuid. In general, in this ID setting, the zero- 
momentum occupation number is a poor diagnostic of su- 
perfiuid transitions, since it may simply refiect a change 
in the value of the critical exponent within a superfiuid 
phase, rather than the onset of exponential correlations. 
Direct evaluation of the correlation functions (a^(a;)a(0)) 
and {rn) {x)rn{fi)) provides a clearer picture in ID, espe- 
cially in the case of a finite size system. Our results are 
fully consistent with the absence of a purely atomic su- 
perfiuid phase in this region of the phase diagram. This 
is compatible with the predictions of mean field theory 
in higher dimensions [52^34j . 

In addition to the zero-momentum occupation num- 
bers, the authors of Ref. [37] also consider the visibil- 
ity. The visibility is related to the momentum occupation 
numbers (22) via [5S] 



^(fc) 



(24) 



where n"^^^ (n""") is the maximum (minimum) in the 
momentum space occupation number distribution. In the 
present context this is identified as 



»a(0) - ng^n) 



(25) 



In a superfiuid phase where ?^q(0) diverges with increas- 
ing system size, the visibility Va approaches unity as 
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FIG. 11. Atomic and molecular correlation functions within 
the AC+MC phase obtained by DMRG with periodic bound- 
ary conditions. We set = —2 and \J — 1.5. (a) Atomic 
correlation function (ajoi+r) as a function of the reduced sep- 
aration r/L for different system sizes, (b) Normalization fac- 
tor Af obtained from panel (a) using Eq. ( |26[ ). (c) Correla- 
tion exponent a obtained from panel (a) using Eq. (261. 



Rescaling the data in panel (a) using the extracted exponent 
a leads to data collapse. This confirms the applicability of 
the conformal result (26 1 within the AC+MC phase. This 
corresponds to power- aw atomic correlations for separations 
r > 3ao. The remaining panels show the corresponding re- 
sults for molecules. 



L — > oo. In Ref. [37 it was argued that the molecu- 
lar visibility within the AC-I-MC phase failed to saturate 
at this value close to the MI boundary. In order to gain 
a quantitative handle on this issue we need to exploit 
the finite-size dependence of the superfluid correlations 
within the AC+MC phase. In a system with periodic 
boundary conditions the two-point function of a primary 
field 0{r) at position r can be obtained by conformal 
transformation 86] : 



Lsin(^) 



(26) 



where a is the critical exponent in the thermodynamic 
limit, r = I?'! — f2| is the separation, and A/" is a con- 
stant pre-factor. It follows that the rescaled combination 
L°-{0{ri)0{r2))L is a prescribed scaling function of the 
reduced separation r/L. The confirmation of this behav- 
ior for the atomic and molecular correlation functions 
within the AC+MC phase is shown in Fig. 11 Given 
this agreement we may substitute the conformal result 
(261 into Eq. (22 1 in order to obtain formal expressions 



for the finite-size dependence of the atomic and molecu- 
lar visibilities in Eq. (24). In Fig. 12 we show the results 



of this conformal extrapolation, where we further incor- 
porate the exact DMRG results for the short distance 
behavior with r < 3ao, where oq is the lattice spacing. It 
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FIG. 12. Finite-size scaling of the atomic and molecular visi- 
bilities within the AC+MC phase. The circles correspond to 
Va and the squares to V,n obtained by Fourier transformation 
of the correlation functions obtained by DMRG. The crosses 
and stars correspond to Fourier transformation of the con- 
formal result ( 26 1 supplemented by exact DMRG results for 
the correlators at small separations r < 3ao. The solid line 
indicates the results of conformal extrapolation (described in 
the text and justified by the scaling collapse in Fig. 11 1 sup- 



plemented by the exact DMRG results for small separations 
r < 3ao. (a) With = —2 and U — 1.5 both Va and Vm 
extrapolate to unity in the thermodynamic limit, (b) Close 
to the MI transition with e,„ = — 2 and U = 2.2 both Va and 
Vm approach unity as L — >■ cxa. This is in direct contrast to 
naive polynomial extrapolation (dashed) which erroneously 
suggests that the molecular visibility is less than unity. 



is readily seen from the solid lines in Fig. [12] that both 
the atomic and molecular visibilities extrapolate to unity 
in the thermodynamic limit. In particular, close to the 
MI boundary there are strong deviations from the results 
that would be obtained by naive polynomial extrapo- 
lation as indicated by the dashed lines. In Fig. [T3|(a) 
we use the conformal extrapolation procedure to track 
the atomic and molecular visibilities within the AC+MC 
phase. The results are consistent with unity right up to 
the MI boundary. For comparison, in Fig. 13 b) we show 
the results that would be inferred using a naive polyno- 
mial extrapolation. The results are in accordance with 
those of Ref. |37| , but differ markedly from the asymp- 
totic visibilities obtained by conformal extrapolation as 
shown in Fig. 



13 



To summarize the results of this section, within the 
AC+MC phase the finite-size dependence of the atomic 
and molecular momentum space diagnostics is in com- 
plete agreement with power-law correlations for both the 
atoms and the molecules. This behavior persists right up 
to the MI boundary and provides further evidence for the 
absence of a purely AC phase. This is analogous to ex- 
pectations in higher dimensions arising from mean-field 
theory analyses 
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FIG. 13. (a) Atomic and molecular visibilities Va (circles) and 
Vm (crosses) within the AC+MC phase for = —2 obtained 
by DMRG with up to L = 64 and periodic boundaries. We use 
the conformal extrapolation procedure described in the text 
in order to obtain the asymptotic results as L — >■ cxa. Both 
Va and Vm are unity right up to the MI boundary, indicating 
the presence of both atomic and molecular superfluidity, (b) 
Naive polynomial extrapolation erroneously suggests that the 
molecular visibility is less than unity in the AC+MC phase. 



VI. ENTANGLEMENT ENTROPY 

Having established good agreement between field the- 
ory and DMRG for the MC and AC+MC phases, let 
us now examine the quantum phase transition between 
them. A key diagnostic in this ID setting is the central 
charge, c, which is a measure of the number of critical 
degrees of freedom. This may be obtained from the en- 
tanglement entropy. For a block of length I in a periodic 
system of length L, the von Neumann entropy is given 
by Sl{1) = -Tri{piliipi), where pi = Ti'L-iip) is the 
reduced density matrix. One obtains |87l 



5l(0= gin 







— sm 









Si 



(27) 



where si is a non-universal constant and where the cor- 
rections are small when the chord length is large fHS?- 
195] . As may be seen in Fig. 14 |a), the numerically ex- 
tracted central charge of the MC phase yields c = 1, as 
one would expect for a single free boson, with coexist- 
ing gapped degrees of freedom; the adjacent panel shows 
the same results plotted against the conformal distance 
I = ln[(L/7r) sin(7r//L)] in order to yield a linear slope of 
c/3. It may be seen from Fig. ^c) that the AC+MC 
phase also has c = 1. Note that it is not c = 2 as would 
be the case for two independent Luttinger liquids. This 
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FIG. 14. Entanglement entropy Sl(1) obtained by DMRG 
with 1/ = 64 and periodic boundaries. We consider horizontal 
scans through Fig. [l] with U = 0.5 and V = 0.7. (a) Within 
the MC phase with = —3.5 we find c « 1 corresponding to 
a gapless superfluid. (b) In the vicinity of the MC to AC+MC 
quantum phase transition we find c « 3/2. This corresponds 
to the presence of additional gapless Ising degrees of freedom 
coexisting with superfluidity, (c) Within the AC+MC phase 
with Em = —4 we find c ~ 1 corresponding to an effective free 
boson. The panels on the right correspond to the same data 
as on the left, but are plotted against the conformal distance 
/ = ln[(L/7r) sin(7rZ/_L)] in order to yield a linear plot with 
slope c/3. The offset between the different curves within each 
panel is due to the non-universal contribution in Eq. (27 1. 



reflects the coupled nature of the atomic and molecu- 
lar condensates in the AC+MC phase, with additional 
gapped Ising degrees of freedom; the Feshbach term is 
relevant and drives the Z2 sector massive. Close to the 
MC to AC+MC transition, where the anticipated Ising 
gap closes, one expects the central charge to increase to 
c — 3/2, due to additional critical Ising degrees of free- 
dom with c = 1/2. This is confirmed by our DMRG 
simulations in Fig. 14 h) . Further support for this Z2 



transition is obtained from the difference 



AS{L) = Sl{L/2) - Sl/2{L/4) = - ln(2) + 



(28) 



as a function of e„j. For a given system size this displays a 
peak, whose location coincides with the MC to AC+MC 
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FIG. 15. Entanglement entropy difference AS{L) on transit- 
ing across the phase boundaries shown in Fig. [T] We use pe- 
riodic boundaries with up to L = 64 and work away from the 
multicritical point, (a) The transition from MC to AC-fMC 
yields c ~ 3/2 corresponding to an Ising transition coexisting 
with a gapless superfiuid. (b) The transition from MC to MI 
yields c « 1 and is consistent with a molecular KT transition, 
(c) The transition from AC-I-MC to MI with = 4 yields 
c « 1 and is consistent with an atomic KT transition, (d) The 
transition from AC-I-MC to MI with em = —3.5 appears to be 
compatible with the approach towards c « 1 with increasing 
L, although the finite-size effects are stronger than those in 
panel (c). Panels (a) and (c) are adapted from Ref. [63] . 



quantum phase transition obtained via the vanishing of 
the single-particle gap, Eig = 0, as shown in Fig. [TJ see 
Fig. |15[ a). The evolution with increasing system size is 
consistent with the passage towards c = 1 in the super- 
fluid phases, and c — S/2 in the vicinity of the transition. 
This behavior may be contrasted with that observed at 
the superfluid-Ml transitions in Fig. [l] away from the 
multicritical point. As may be seen in Fig. |15[ b), in 
passing from the MC phase to the Ml, the central charge 
remains pinned at unity. This is consistent with a KT 



transition for the molecules. Likewise, in passing from 
the AC-I-MC phase to the Ml, we find c = 1 again; see 
Fig. 15 ^c). This is consistent with a KT transition for 
the atoms. We have checked that this atomic KT be- 
havior persists in Fig. [l] up to a value of e„i — —3.5; 
see Fi g. [l5[ d). It is notable that the finite-size effects 



15 



in Fig. |15[ d) are much stronger than those in Fig. 
although both are compatible with c 1 at the Ml tran- 
sition. A detailed analysis of the multicritical region in 
Fig. [l] requires further investigation. 



VII. ISING SCALING REGIME 

Having provided evidence for a Z2 quantum phase 
transition occurring between the MC and AC-j-MC su- 
perfiuids, we now demonstrate how to extract both the 
Ising order parameter, (</>), and the Ising correlation 
length, ^, in the presence of the additional superfluid 
degrees of freedom with c = 1. 



A. Ising Correlation Length 

The Ising correlation length, ^, may be obtained 
from the atomic and molecular correlation functions dis- 
cussed in Sec. |lVBl Within the Z2 disordered MC phase 
the atomic correlations {a'^{x)a{0)) ~ a;~''™/''Ko (x/f ) 
decay exponentially, whilst the molecular correlations 
{m^{x)m{0)) ~ a:""'" decay as a power-law. At a 
given point in parameter space we may use the molec- 
ular Green's function to determine the exponent Vrm 
and thereby extract the Ising correlation length from the 
atomic correlations. This approach is outlined in Fig.flG] 
In the vicinity of an Ising quantum phase transition one 
expects that ^ \M — A4c\'^ where v = 1 is the Ising 
correlation length exponent and Ai is a suitable mass 
scale parameterizing the departure from criticality. Un- 
fortunately, it is non-trivial to express A4 in terms of the 
microscopic parameters of the lattice model ([T]) . A naive 
analysis gives ^A ^ ko + Kipm + ^2 y/Pm , where pm is the 
density of molecules, and kq ~ Eq, ki ~ Uam, ^2 '■^ 2g 
are constants. Expanding the square root according to 
\/P^ ~ y/f4n + (Prn ~ Pm)/VP?n suggests that Sufficiently 
close to the Ising transition 

r'^lPrn'P'J. (29) 

This Ising behavior with 1/ = 1 is confirmed in Fig. |16| 



B. Ising Order Parameter 

In the ordered phase of the Ising model ^ we have 
{(j)) ^ \A4 — -Mcl^ where (3 = 1/8 is the Ising magneti- 
zation critical exponent. From the discussion above one 
thus expects that 
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FIG. 16. DMRG results used to extract the Ising correlation 
length ^ within the Z2 disordered MC phase with L — 256 and 
open boundaries, (a) With em ~ —3.85 the molecular corre- 



lation function {mlmj) ~ 0.81 1 



decays as a power- 



law, (b) Using the previous exponent Um ~ 0.18 we extract 
the Ising correlation length ^ « 43.6 from the exponential de- 
cay of {a\aj) ~ |z — j|~''™/*Ko(|i — (c) Repeating the 
above procedure we plot (circles) versus the departure 
of the molecular density pm from its value at the MC to 
AC-I-MC transition. Close to the transition the results are in 
good agreement with the Ising relation ~ |Pm — /OmT with 
1^ = 1. The triangles correspond to extracting ^ directly from 
the ratio n{\i - j\) = {a\a,)^ / {mlm,) ~ [Ko(li - 



where Pm is the density of molecules. In order to test 
the validity of Eq. ( 30 ) , we must first extract the Ising 



order parameter from a finite-size scaling analysis of the 
atomic correlations. As foUows from Eq. (16), within the 
Z2 ordered AC-I-MC phase one has 



(31) 



where ,4 is a normalization amplitude. In Fig. |17[ a) we 
show DMRG results for (a^(x)a(O)) in the vicinity of the 
MC to AC+MC quantum phase transition. A direct fit 
to Eq. (31) yields A{4>)'^ . Repeating this procedure for 
different system sizes one obtains an estimate for A{(t>)'^ 
in the thermodynamic limit; see Fig. [iTjjb). In Fig. [Is] 
we show the variation of this order parameter with the 
molecular density. The results are in good agreement 
with the theoretical prediction in Eq. ( 30 ) and the Ising 
critical exponent (3 = 1/8. 



C. Correlation Function Ratio 

In the above discussion we have extracted the Ising cor- 
relation length and the Ising order parameter through a 
direct finite-size scaling analysis of the atomic and molec- 
ular correlation functions. An alternative approach is to 
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FIG. 17. (a) Atomic correlation functions within the AC-I-MC 
phase obtained by DMRG with up to L = 512 and open 
boundaries. We set U = 0.7 and consider e m ~ —3.75 (left) 
and em = -3.79 (right). A direct fit to Eq. (3l]) yields A{(t>)^ 
for each value of em, where {(jj) is the Ising order parameter 
and ^ is a non-universal constant prefactor. Changing the 
fitting interval gives an estimate of the error bars, (b) Ex- 
trapolation of A{4>)'^ to the thermodynamic limit using linear 
extrapolation of the largest three system sizes is indicated by 
the solid line. An estimate of the error bars in the thermo- 
dynamic limit is obtained by comparing to a quadratic fit of 
the data shown by the dashed line. These results are plotted 
as a function of the molecular density in Fig. [18] in order to 
confirm Ising behavior with /3 = 1/8. 



consider the behavior of the ratio 

_ (at(x)a(0))4 



TZ{x) 



(mt(x)m(G)) 



(32) 



in analogy to the considerations of Refs. [HI [33] for the 
confinement-deconfinement transition of Cooper pairs in 
ID fermion systems. 

In the Z2 disordered MC phase {a^{x)a{0)) ^ 
x~'''^/^Ko{x/£_) and (TO^(a;)m(0)) - x-"--. It follows 
that the power-law prefactors cancel out in this ratio: 



nix) ^ [Ko(a;/e)]' 



(33) 



As such, this ratio should exhibit exponential decay in 
the MC phase. This is confirmed by our DMRG results 
in Fig. ^a). A direct fit to Eq. ([SSf yields the Ising 
correlation length shown by the triangles in Fig. |16[ 
In the Z2 ordered AC-I-MC phase {m^x)m{0)) 



and (al'(x)a(O)) 
given by Eq. (fl7| 



((/)(a;)0(O))x-'''"/4 where ((/)(a;)0(O)) is 
It follows that 



nix)^{(b)'[l + n-'F{x/Oy 



(34) 
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FIG. 18. DMRG results for the Ising order parameter in the 
Z2 ordered AC+MC phase with up to L = 512 and U = 0.7. 
(a) Variation of the Ising order parameter squared A{4')'^ 
where ^ is a constant prefactor, versus the deviation of the 
molecular density pm from its value at the MC to AC+MC 
quantum phase transition, (b) Variation of A'^{(j))^ versus the 
molecular density difference. The results are in good agree- 
ment with the Ising magnetization relation (0) ~ \pm — Pm\^^^ 
with P = 1/8. The inset shows analogous results obtained 
from the plateau value of TZix) for L = 512, as indicated in 
Fig. 19 c). The error bars are estimated from the magnitude 
of |7l(s = 128, L = 512)-7^(a; = 64, L = 256) |. (c) Finite-size 
scaling of the thermodynamic molecular density used in pan- 
els (a) and (b). Panels (a) and (b) are adapted from Ref. [63] . 



where F{z) is given by Eq. ( 18 ). At leading order TZ{x) ^ 
{(j))^ and one thus expects TZ{x) to develop a constant 
plateau that is proportional to the Ising order parameter. 
This is confirmed by our DMRG results in Fig. 19 c). 

In addition to these results for Tl{x) which are valid 
within the superfiuid phases, one may also explore the 
vicinity of the quantum phase transition between them. 
At the Ising critical point (0) = but 



Wx)m) - (- 

V X 



(35) 
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FIG. 19. Ratio 7^(j^ - j\) = (alaj)"^ /{mlmj) of the atomic 
and molecular correlation functions with U = 0.7. (a) In the 
Z2 disordered MC phase with tm = —3.85 the ratio 7l{\i — 
j\) ~ [Ko(!j — jl/O]"^ exhibits exponential decay, (b) In the 
vicinity of the MC to AC-I-MC quantum phase transition the 
ratio TZ{\i — j\) ~ 1/1* — i I decays with a universal power- law 
corresponding to the Ising critical exponent rj = 1/4. The 
line is a fit to 7?. = Ao\i — over the interval S < \i — j\ < 
48 for L = 512 with Ao ~ 0.011 and Ai « -1.01. (c) In 
the Z2 ordered AC-I-MC phase with = —3.75 the ratio 
^(N ~ il) ~ {4')^ exhibits a plateau corresponding to a non- 
zero Ising order parameter. 



decays as a power-law where 77 = 1/4 is the Ising pair 
correlation exponent. It follows from Eqs. (|6| and ^ 



that the atomic Green's function at criticality is given by 
(at(x)a(O)) 



(mm) 

ao 

X 



.t3(k) -13(0) 



(36) 



On passing from the Z2 ordered AC-I-MC phase towards 
the Ising quantum phase transition, the power-law decay 
of the atomic Green's function is therefore enhanced by 
77 = 1/4 due to the presence of additional critical Ising 
degrees of freedom. It follows that 



n{x) 



Oq 
X 



(37) 



exhibits universal power-law decay in the vicinity of the 
MC to AC-fMC quantum phase transition. This is con- 
firmed by our DMRG results shown in Fig. 19 ^b). This 
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provides direct evidence for the Ising correlation expo- 
nent, ?7 = 1/4. These resuhs demonstrate that the ratio 
TZ{x) may be used to explore both the critical and ofF- 
critical Ising behavior at the MC to AC+MC transition. 
The characteristic signatures of TZ{x) parallel those ob- 
served in Refs. [42l|43] for the confinement-deconfincmcnt 
transition of Cooper pairs in ID fermion systems. 



VIII. CONCLUSIONS 



In this manuscript we have explored the phase diagram 
of bosons interacting via Feshbach resonant pairing in a 
ID optical lattice. We have presented a wide variety 
of evidence in favor of an Ising quantum phase transi- 
tion separating distinct paired superfluids. We have also 
provided a detailed characterization of these phases, in- 
cluding the behavior close to the Mott insulating phase 
boundary. For the investigated parameters, our DMRG 
results are consistent with an Ising quantum phase transi- 
tion approaching both a molecular KT transition and an 
atomic KT transition. This is compatible with mean field 
theory predictions for the continuum model in higher di- 
mensions. However, recent results for pairing phases in 
a 2D classical XY model suggest the possibility that the 
Ising transition may over-extend beyond the multicritical 
point [53]. In view of this possibility, in a distinct but 
closely related model, it would be profitable to explore 
the multicritical region in more detail. A clear verdict on 
this issue for the present ID quantum model requires fur- 
ther analytical and numerical investigation and we will 
return to this question in future work. It would also be 
interesting to explore the phase diagram for a broader 
range of parameters, with a specific focus on the choice 
of atomic species and experimental constraints. Even in 
the presence of strong three body losses, the emergent 
phase diagram may exhibit notable similarities [451151) . 
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Appendix A: Hilbert Space Truncation 

Throughout the main text we truncate the local 

Hilbert space to allow up to a maximum of ti^ = 5 atoms 

and rirn, — 5 molecules per site. In the regime of large 
-4.46 1 ' 1 ' 1 ' 1 ' 
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FIG. 20. DMRG results for a horizontal scan through Fig. [l] 
with 1/U — 2.5. We consider up to L = 128 and extrapolate 
to the thermodynamic limit, (a) Evolution of the ground 
state energy Eq with increasing local Hilbert space restriction, 
(b) Evolution of the excitation gap Eig with 
increasing rir, showing very little change beyond = 5. 



t/U, where inter-site hopping is strongly favored, one 
should check the validity of this approximation. Here we 
discuss the evolution of physical observables with increas- 
ing Hilbert space restriction parameter, rir — Ua = rim- 
For the largest value of l/U = 2.5 used in Fig. [T] the 
results converge with increasing n^. For example, in 
Fig. [20|^a) we show the evolution of the ground state en- 
ergy Eq with increasing n,.. The results show very little 
variation beyond Ur = 5. Likewise, in Fig. 20 b) we 
monitor the excitation gap Eig with increasing Ur. The 
data again show very little change beyond rir = 5. The 
associated MC to AC-I-MC phase boundary in Fig. [l] is 
therefore robust to increasing n^. In a similar fashion, in 
Figs. [21] and [22} we examine the evolution of the atomic 
and molecular correlation functions. The results show 
clear convergence in both the MC and AC-l-MC phases. 
The excellent agreement between our DMRG results and 
field theory predictions also lends a postiori support for 
this level of Hilbert space restriction with rt^ = 5. 
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FIG. 21. (color online). DMRG results with L = 128 and 
U — 0.5 within the MC phase shown in Fig. [T] We show the 
evolution of the atomic and molecular correlation functions 
with increasing local Hilbert space restriction n,.. We set 
e,n — —4 {em = —6) in the upper (lower) panels. 




FIG. 22. (color online). DMRG results with L = 128 and 
U = 0.5 within the AC+MC phase shown in Fig. [l] We 
show the evolution of the atomic and molecular correlation 
functions with increasing local Hilbert space restriction n^. 
We set em ~ —2 {em ~ —3) in the upper (lower) panels. 
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